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Abstract 

We study the final distribution of the winding numbers in a ID superconducting ring that is 
quenched through its critical temperature in the absence of magnetic flux. The study is conducted 
using the stochastic time-dependent Ginzburg-Landau model, and the results are compared with 
the Kibble-Zurek mechanism (KZM). The assumptions of KZM are formulated and checked as 
three separate postulates. We find a characteristic length and a characteristic time of this system. 
Besides the case of uniform rings, we examined the case of rings with several weak links. For 
temperatures close or below Tc, the coherence length does not characterize the correlation length. 
In order to regard the winding number as a conserved quantity, it is necessary to allow for a short 
lapse of time during which unstable configurations decay. We found criteria for the validity of 
the ID treatment. The is no lower bound for final temperatures that permit ID treatment. The 
variance of the winding number, (ra^), depends not only on the quenching time tq, but also on the 
size of the change of the reduced temperature. For moderate quenching times, we find the scaling 
{v?) oc "^Q^^^j as predicted by KZM in the case of mean field models; if tq is smaller than the 
typical time over which unstable configurations decay, the dependence is weaker. 

PACS numbers: 74.40. +k, 74.81.Fa, 05.70.Fh, ll.15.Ex 
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I. INTRODUCTION 



The Kibble-Zurek mechanism^"- (KZM) aims at the description of kinetically governed 
phase transitions induced by rapid decrease of the temperature. KZM may be didactically 
expressed as the division of the transition into three stages, with underlying postulates for 
each stage. During the initial stage the system is in a high symmetry phase; it explores lower 
symmetry states due to thermal fluctuations, but is symmetric on the average. The system 
is able to follow the change of the temperature and may be considered to remain in thermal 
equilibrium. During the intermediate stage the critical temperature is crossed. During this 
stage the system is sluggish, so that its final state will be very similar to the initial state. 
In the final stage several lower symmetry states are long lived; the selected state will have 
symmetry properties that are determined by the state that was randomly occupied by the 
system when it entered the second stage. 

KZM has been tested by many experiments and simulations, and various degrees of agree- 
ment have been obtained. We will be interested in the case of systems with loop topology.-""— 
In this paper we consider a superconducting loop cooled at a fast rate from slightly above 
to slightly below its critical temperature, in the absence of applied magnetic field, and 
investigate the possibilities for spontaneous emergence of a permanent supercurrent. 

II. QUANTITATIVE FORMULATION OF KZM POSTULATES 

The first postulate specifies the temperature at which the system becomes sluggish: 

Postulate 1 The crossover from the initial to the intermediate stage occurs when r(e) = 
|e/e|. Here e = {Tc — T)/Tc is the reduced distance from the critical temperature, e is the 
derivative of e with respect to time, and r(e) is the relaxation time of the system. 

The following postulates are formulated here for the case of a system with loop geometry, 
described by a complex order parameter. The moment at which e = will be taken as t = 
and the moment at which the system passes from the initial to the intermediate stage will 
be denoted as t = —i. 

Postulate 2 We denote by ^ the coherence length of the order parameter at t = —t and 
by C the perimeter of the loop. The behavior of the system can be estimated by enviewing 
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it as if it were divided into M = C pieces, such that in each piece the phase of the order 
parameter is uniform, and such that there is no correlation between the phases in different 
pieces. 

Let us denote by 9i the phase difference between regions i and i + 1. It follows from 
Postulate [2] that the variance of this phase difference is {Of) = ttVS. Let Oc = T.'ti be 
the phase difference accumulated around the loop. Since {OiOj) = for i ^ j, if we ignore 
the constraint that 6c has to be an integer multiple of 2tt, we obtain 

{el) = ATvrVS . (1) 

It turns out that the constraint has no effect. Taking the ^j's uncorrelated for i = 1, . . . , Af—1 
and picking 6'^ as the phase difference with the smallest absolute value such that n = dc/27r 
is integer, Eq. ([1]) is recovered for A/" > 3. n is the winding number of the order parameter, 
and in view of Eq. ([T]) obeys 

(n^) = Af/12 . (2) 

Postulate 3 The winding number can be regarded as a topological invariant that remains 
unchanged for t > —t. 

This postulate is far from plausible. Although n is a discrete variable, there is no impediment 
for a continuous passage of the order parameter to a different winding number, provided that 
it vanishes at some point. 

The idea behind Postulate 3 is that during the period —i<t<i the system is sluggish 
and its order parameter practically does not change, and after t = t there is an energy 
barrier between states with different winding number. This statement would be true if at 
t = t the order parameter were close to local equilibrium, which most probably will not be 
the case. Since the order parameter is random and nearly vanishes at t = —i, the same is 
true at t = t, and we have no good reason to expect a considerable energy barrier for the 
passage to a different winding number. 

Postulate 3 has been amended in later publications (e.g. Ref. ilOj) . The revised claim is 
that (i) there is a time t^ at which the order parameter shoots up sharply; (ii) n remains 
constant after t^; (hi) tAOci and tA is of the order of i; (iv) since the coherence lengths are 
the same for e and — e, it follows from (i)-(iii) that the scaling properties of the system will 
be the same as if Postulate 3 were true. 
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Most of the literature on KZM studies the scahng of the density of defects (the winding 
number in our case) with the rate at which the system is cooled. For this purpose, Postulate 
1 is of central importance. In the present study we will first consider the case in which the 
system is cooled at an ideally fast rate and will check the applicability of Postulates 2 and 
3; verification of Postulate 1 will be postponed to Sec. |Vll 



III. OUR MODEL AND SYSTEM 



A simple approach for the description of the dynamics of a superconducting sample is the 
time-dependent Ginzburg-Landau model, with the addition of Langevin terms that bring 
thermal fluctuations into account (e.g. Ref. Illl ): in the case of ID systems, simplifications 
are possible. 12 We assume that the magnetic field induced by the current around the loop 
is negligible and write the Ginzburg-Landau energy as 
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where m is the mass of an electron pair, s is the arc length and the integral covers the 
loop, ^(0) is the coherence length at T = 0, is the order parameter, A is the tangential 
component of the electromagnetic vector potential, $o is the quantum of flux, w is the 
cross section of the loop and /3' = 'ine^n? /mc^, where e is the electron charge, k is the 
Ginzburg-Landau parameter and c the speed of light. The usual coherence length in the 
Ginzburg-Landau model is ^(e) = ^(0)/|e|^/^. We assume that i}) and A only depend on 
time and s, and not on the lateral position. A is required in the energy functional in order 
to take account of the fluctuations of the electric field. 

In the numerical procedure, the loop is divided into N segments and the integral is 
approximated by the sum^^ 
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where the segment + 1 is identified with the segment 1, and the evolution of ipj and Aj 
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during a period of time At ^ r(e) is given by 



dG 
dG 

A, (t + At) = A, (t) - Ta,j ^At + r]A. (5) 
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Here F^j = NmDj/nrGwj, = Nc /ajCwj, D is the diffusion coefficient and a is the 
conductivity, r^i, 772 and rj^ are random numbers with gaussian distribution, zero average 
and variances {rjl) = {rj^) = 2r^jkBTAt, (77^) = 2T AjksT At, where ks is the Boltzmann 
constant. 

The material parameters in the model can be evaluated in terms of Tc, the Fermi wavevec- 
tor kp, the mean free path le and the electron density Ue- Using BCS, dirty limit and free 
electron gas approximations. 

An important characteristic length is^ = (w^l/32TT^K'^kBTcY^^ . For e ^ 0, is 
the order of the length over which ip is not expected to vanish for a fluctuation such that 
G ~ ksT,. 

In order to have a ID situation, and ^(e) should be larger than the linewith and the 
thickness. In the absence of significant magnetic fields, the magnetic penetration depth 
n^{e) is unimportant. 

We note that for the present model the relaxation time does not diverge at T^. First, 
relaxation of the electromagnetic potential is insensitive to the distance to T^, and indirectly 
drives the phase of ip- Moreover, for e = the remaining terms in Eq. ([3]) lead to a residual 
relaxation time of the order of ^"i/D- 



IV. VERIFICATION OF POSTULATE 2 



Since this postulate is intended to be applied at a moment such that the system is still 
able to follow the temperature variation, in this section we consider systems that are allowed 
to reach thermal equilibrium at some fixed reduced temperature e = (Tc — T) /Tc, and check 
the result predicted by Eq. ([2]). 

We studied rings of length G = 0.1cm and electron density Ue = lO^^cm"^, divided 
into 300 computational cells. The system evolved during a lapse of time 2^'^/D, in steps of 
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FIG. 1: Variance of the winding number of the order parameter for a uniform ring close to T^. 
For T > Tc (resp. T < Tc) the abscissa is C/^ (resp. — C/^), i.e., the abscissa is an increasing 
function of the temperature. The short horizontal lines that cut the y-axis mark the values 
□ : w = 5 X lO-'^cm^, 4 = IQ-^cm, = IK; *: w = lO^'^cm^ 4 = lO^^cm, = 5K; o: 



w 



IQ-'^cm^, 



2.5 X 10-'^ cm, Tc = 5K;o: w = 2.5 x IQ-^cm^, 



lO-^cm, Tr = 20 K. 



3.3 X 10^^^|/-D. At the end of this evolution the winding number n = (l/27r) ^ a.Tg{ip*il!j^i) 
was evaluated. For each temperature above (respectively below) Tc this procedure was 
repeated 800 (respectively 400) times. The evolution time seemed to be sufficient, since 
there was no significant correlation beween the initial and the final winding number; besides, 
except for temperatures considerably below Tc and unrealistic initial distributions, longer 
runs lead to essentially the same results. 



A. Uniform loop 

We studied both the case T > T^. and the case T < T^.. The initial distribution of the 
order parameter was taken from a quadratic approximation, as described in the Appendix. 
Figure [T] shows our results. We observe that (i) {n^) is a universal function of ^ and ^p] the 
rhombs and the circles in the figure, that were obtained for samples with different parameters 
but the same value of lie along the same curve; (ii) for T ^ Tc, {n^) is surprisingly close 
to C/12^^; (iii) for T > Tc and ^ ^ {n^) does not depend significantly on ^^; (iv) as 
the temperature decreases below Tc, (n^) decreases until the probability for n 7^ becomes 
negligible. 
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FIG. 2: Variance of the winding number as a function of the ratio between the length of the loop 
and the correlation length. Insets: correlation length as a function of e = 1 — T/Tc- Parameters 
and symbols as in Fig. [TJ 

The fact that below Tc {n^) decreases as C/^{e) increases indicates that ^(e) is not 
necessarily the characteristic length over which the phase remains essentially unchanged. 
The relevant length ought to be the correlation length, which can be defined as follows: 
we define the autocorrelation function K{s) = {Re[ip*{s')ip{s' + and then the 

correlation length sk is obtained from a fit K{s) ~ exp(— s/s^-)- 

sk and ^(e) are not equivalent concepts for two reasons. First, .^(e) arises from a compe- 
tition between the stiffness and the condensation energy, whereas in the case of sk also the 
thermal energy and the quartic term enter the competition. More fundamentally, ^(e) is a 
healing length, which is relevant when ip has to obey constraints, e.g., it has to vanish at 
every vortex; in the problem we are considering, the only constraint is periodicity. 

The lower inset in Fig. [2] shows the temperature dependence of sx(e), scaled by ,^/3 (for 
two samples), and the higher inset compares the temperature dependences of Sxi^) and ^(e). 
Sufficiently above sk ~ Near Tc, sk{^) does not diverge; = 0) is of the order 
of ^p. Contrary to ^(e), sx(e) is a monotonic function and rises sharply as the temperature 
decreases below Tc. If sk is not significantly smaller than C, then the autocorrelation is 
influenced by the connectivity of the loop and K{s) is not well fltted by an exponential 
function. 

The main panel in Fig. [2] shows (n^) as a function of C/sx- Within our numerical 
accuracy we obtain a universal line for both samples, which is fltted, both below and above 
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Tc, by (n^) = 0.106C/sk ~ C/{12 x O.Ssk), i.e., O.Ssk is the effective length over which the 
phase is essentially uniform. 

B. Chain of Junctions 

We now study a loop with Afw weak links. For a situation in which the phase is fairly 
uniform between weak links but uncorrelated across them, we could argue that the present 
case is equivalent to that of Eq. ([2]), with Af = Afw 

We consider a temperature that is above and very close to it. Between weak links we 
took w = 5 X lO^^cm^, = 2 x 10~^cm, Tc = 1 K, e = —2 x 10^''; the other parameters 
were taken as in the previous section. The weak links had lengths iw, had a central part 
with e = —3 X 10^^ (i.e., had lower Tc), and had shorter mean free paths. The outer parts 
of the weak links had positive values of e, chosen to compensate the influence that the weak 
links would have on the rest of the loop due to proximity. This time we took iljj{t = 0) = 0. 
In the evaluation of the winding number, the weak links were skipped. 

The dashed line in Fig. [3] corresponds to the ideal case Af = Afw Taking into account that 
the weak links are not perfect would lead to = pAfw, where p G (0, 1) is some "opacity 
coefficient." We also have to take into account that, according to the previous section, along 
the length C — Afw^w of the strong part of the loop, (C — Afwiw)/^i3 divisions into different 
regions effectively occur. Therefore, instead of Eq. ([2]) we expect 

{n') = C/12^p + (p - C/e^)AC/12 . (7) 

The moderately weak links in Fig. [3] had length = 1.67 x 10^'^ cm, mean free path 
10~^cm, and the region with low Tc had length 3.33 x 10~^cm; for these weak links we 
obtained p = 0.45. The weaker links in the figure had length 3 x 10~'^cm, mean free path 
4 X lO^^cm, and the region with low Tc had length 10~^cm; for these weak links we obtained 
p = 0.96. 

V. VERIFICATION OF POSTULATE 3 

The equilibrium distributions that were obtained in Sec. IIVI will now be taken as initial 
distributions, and we will study under what conditions the winding number remains un- 
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FIG. 3: Variance of the winding number for a loop with weak links. The dashed line corresponds 
to Af = Myj. •: moderately weak links; □: very weak links. The solid lines corresponds to Eq. ([7|). 



changed if the system is instantaneously cooled to a new temperature Ty = Tc(l — tf) and 
left at this temperature for a considerable time. 



A. Pure ID case 



Since ^(e/) may be shorter than the linewidth or the thickness of the loop, our first 
concern should be whether the system can still be treated as ID. Following our findings in 
Sec. IIV Al we may expect that the relevant requirement is sxiej) ^ w^^'^. We will assume 
here that the ID treatment is still applicable; deviations from ID will be studied in the 
following section. 

It is known that in the absence of magnetic flux a state with winding number n such that 
> Cy2n^e{^) + 1 is locally unstable (e.g. Ref. Q), so that Postulate 3 necessarily fails 
in this situation. In view of Postulate 2 and Eq. ([2]), this situation can be avoided provided 
that e^(e/) < Ci. 

We should require that no further changes of the winding number occur after a state 
close to a local minimum of G is reached. For this effect, the energy barrier for the 
passage to another winding number should be significantly larger than ksTf. It fol- 
lows from expression ([3]) that, near a local minimum, the energy per unit length is 
-(/i2y;/8m/3')(l/^^(e/) - 47rV/C2)2 ^ -h'^w /8m(3'^'^{ef). In order to change the wind- 
ing number, the order parameter has to vanish at some point and the smallest possible 
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energy increment involved is obtained when the order parameter becomes small in a region 
of length ~ ^(e/)- From here, h'^w/8'mP'C,^{ef) has to be significantly larger than ksTj, a 
condition equivalent to the requirement that ^(e/) be significanly less than 

Provided that at time (i) the sample is described by an equilibrium distribution with 
temperature Tc(l — e), (ii) its winding number is n 7^ and (iii) the temperature is instan- 
taneously lowered to Tc(l — e/), we define the "conservation probability" V by 

V{e,ef,t) = Probability that the winding number remains invariant 

during the lapse of time between and t . (8) 

The condition of invariance is required in a strict sense, i.e., n{t') = n{0) for all < t' < t. 
According to Postulate 3, V ought to be 1. 

The markers in Fig. H] show the conservation probability for the sample represented by 
squares in Fig. [H as a function of the initial reduced temperature e. In the case of the black 
markers, ^(e/) ~ 0.02,^/3, so that the final energy barriers are expected to be considerably 
larger than ksTf. We note that V is substantially less than 1. It also appears that, as the 
initial temperature increases above Tc, V decreases exponentially down to practically zero 
within a very narrow range. The red markers are for higher Tf [^{^f) ~ O.OG.^/?]. It would 
seem that the energy barriers are still safely high, but the conservation probability for the 
case €f = 10^^ is quite lower than in the case = 10~^. In order to estimate the amount 
by which V decreases due to thermal jumps over the energy barriers, we repeated some 
simulations with thermal fiuctuations turned off at t = (blue and purple markers); the 
absence of fluctuations leads to conservation probabilities that are larger by about ~ 0.1. 
We therefore conclude that thermal fluctuations are not the main reason for the change of 
the winding number; n usually changes just because the initial order parameter is not close 
to a local minimum of G, and is not within the basin of attraction of the local minimum 
that has winding number equal to n{0). The lower conservation probabilities for small e/ 
could be due to insufficiently fast growth of the order parameter. 

Figure [5] shows the dependence oi V on w (blue and green markers) and on Tc (red and 
black markers), for a range somewhat larger than an order of magnitude that includes the 
sample studied in Fig. HJ with initial temperature close to T^. As a general trend, we may 
regard V as independent of and proportional to w^^^. Extrapolating this result we would 
obtain V{e ~ 0, 10"^, lO^j/D) ~ 1 for w ~ lO^^cm^. Close to the parameters in Fig. HI V 
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FIG. 4: (Color online) Conservation probability V{e,ef, 10(,'^/D) [defined in Eq. dS])] as a function 
of the initial temperature. The red markers are for £/ = 10^'^ and the black markers for e/ = 10~^. 
The blue and the purple markers are obtained if no thermal fluctuations are applied at the final 
temperature. The error bars were estimated assuming a relative uncertainty equal to 1 over the 



square root of the number of counts, it; = 5 x 10 cm 



lO-^cm, Tc = IK. 



is also roughly proportional to C~^/^ and to (not shown). In addition to the general 
trends, we have also noticed (but not systematically investigated) a non monotonic detailed 
structure in the dependence of P on e and on the sample parameters. This dependence may 
be due to the availability of additional channels when the energy landscape varies. 



1. A weaker version of Postulate 3 

At the moment that the sample is (instantaneously) quenched, the order parameter may 
or may not be within the basin of attraction of the local minimum of the energy with the 
same winding number. If it is (respectively, is not), we may regard this situation as a sort 
of "metastable" (respectively "unstable") state. We may expect that unstable states decay 
very fast, so that after waiting a short time the winding number could indeed be conserved. 
With these heuristic ideas in mind, we define [under the same conditions as Eq. ([8])] the 
"restricted conservation probability" Vr by 

Vnle, ef,tu,t) = Probability that the winding number remains invariant 

during the lapse of time between tu and t . (9) 
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FIG. 5: (Color online) Conservation probability as a function of the cross section and of the critical 
temperature of the wire. In all cases C = 0.1cm, if, = 10~^cm and ej = 10~^. When studying 
the dependence on the cross section (lower x-axis), Tc = IK; in this study the blue markers are 
for e = —1.5 X 10~^ and the green markers for e = 1.5 x 10^^. When studying the dependence on 
the critical temperature (upper x-axis), w = 5 x 10~''cm^; in this study the red markers are for 
e = —1.5 X 10~^ and the black markers for e = 1.5 x 10^^. 

A weak version of Postulate 3 (WPS) is then: for an appropriate short time tu, 
VR{e,ef,t,,t)^l. 

Figure E] shows the restricted conservation probability as a function of e for several samples 
and choices for t„. The initial temperatures were taken above Tc. The black markers are for 
tu = 10~^^|/D and the colored markers for shorter waiting times. The squares are for the 
same sample as the squares in Fig. [H all the other markers are for samples in which one of 
the parameters was changed by a factor 10. Figure [6] indicates that WPS is more pervasive 
than could be expected from a non quantitative statement; except for the sample with small 
cross section, VR^e, 10^^, 10^^^|/D, 10^|/D) is indistinguishable from 1 in the entire studied 
range. Note that tu for the black markers is smaller by three orders of magnitude than the 
following period of time during which the winding number is monitored for possible changes 
at the final temperature. 
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FIG. 6: (Color online) Pij(e, 10~^, t^, lO^^/Z?) [defined in Eq. @] as a function of e, for several 
samples and several values of t^. The blue markers are for tu = W~^^'^/ D, the red markers are for 
t„ = IQ-^^j/D and the black markers for tu = IQ-'^^j/D. □: w = 5x lO^^cm^, 4 = lO^^cm, = 
IK, C = 0.1cm; o: w = 5 x lO^^cm^, 4 = lO^^cm, = WK, C = 0.1cm; o: w = 5x lO^'^cm^, 
4 = 10-^cm, Tc = IK, C = 1cm; x: w = 5x lO^^cm^, 4 = lO'^cm, = IK, C = 0.1cm; •: 
w = 5 X 10-^cm^, ie = 10~^cm, Tc = IK, C = 0.1cm. Uncertainty bars are not shown. 

B. Quasi-2D treatment 

1. Criterion for one- dimensionality 



Let us model the cross section of the loop as a rectangle of length 6 (i.e., 6 is the width of 
the loop). How small 6 has to be in order to justify the ID treatment? It is usually assumed 
that 6 has to be smaller than the coherence length and the Pearl length. These would be 
the appropriate requirements if there were a significant magnetic field and the temperature 
were sufficiently far from T^; we claim that in the present problem the relevant requirement 
is 5 < sk, where sk is the correlation length defined in Sec. IIV A[ 

It might be argued that sk was evaluated under the assumption that the system is ID, 
so that adopting the same correlation length for the lateral direction is a sort of circular 
reasoning. As a test for our claim, we consider a simple model in which the lateral direction 
is taken into account. We divide the loop into two halves, each having length C and width 
6/2. We denote the order parameters in the respective halves by ip(i){s) and 4'(2){s) and 
approximate the lateral component of the gradient of ip by 2[ip(^2){s) — 4'(i){s)]/6. Ignoring 
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FIG. 7: (Color online) Lateral correlation of the phase of the order parameter, o; § = sk', *'■ ^ = 



Sample parameters: C = 0.1cm, w = 2.5 x 10 °cm 



10-^cm, Tc = 20 K. 



the lateral component of the electromagnetic potential, which amounts to ignoring Johnson 
noise in the lateral direction, the energy becomes 
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-^(^a)^(2) + ^(2)V'(i)) \wds. 



(10) 



Let us write ipiu) = \4'iu)\^^'p{iXiu)) and consider the average (cos[x(i)(s) — X(2)(s)]) (sta- 
tistical average, for arbitrary s). If (cos[x(i)(s) — X(2)(s)]) ~ 1, it means that the phase 
variation is practically the same in each of the halves of the loop and the ID treatment is 
justified; if (cos[x(i)(s) — X{2){s)]) ~ 0, there is no correlation between the phase variations 
in each of the halves, and we have a qualitatively different problem. We note that the energy 
G'q2d has the same structure as that of a two-band superconductor, and a procedure for the 
evaluation of (cos[x(i)(s) — Xi2){s)]) is available. 

Figure [7] shows the temperature dependence of (cos[x(i)(s) — X(2) (■§)]) for the sample 
represented by circles in Fig. [H When we take 6 = C, (filled symbols), there are temperature 
ranges where (cos[x{i)(s) — X(2){s)]) ~ 0, indicating that a loop of width ^ can not be treated 
as ID; there are also ranges where (cos[x{i)(s) —Xi2){s)]) ~ 1, indicating that wider samples 
could still be treated as ID. On the other hand, taking 6 = sk (empty red symbols) leads to 
intermediate correlations for all the considered cases, showing that this is the verge of the 
ID behavior. 
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There is a moderate tendency for a decrease of (cos[x(i)(s) — X(2){s)]) for 5 = as 
the temperature is lowered, but, in view of our resuhs in the lower inset in Fig. [2|, we 
still conclude that the admissible width for ID treatment does not decrease as temperature 
decreases below T^. It follows that if the order parameter is initially ID, thermalization at 
Tf will not take the system away from one dimensionality, even if ^(7/) ^ 6. 



2. Influence of the induced flux 

In the previous sections we entirely neglected the magnetic field induced by the current 
I{t) that circulates around the loop. For a quasi ID loop, the influence of this magnetic field 
can be taken into account by means of a self inductance L that gives rise to a magnetic flux 
Lcl enclosed by the loop. Self inductance restricts the rate at which current can vary in a 
circuit to the characteristic time = Law/C. Since usually a change in winding number is 
accompanied by a change in the current, our results in Fig. [6] lead us to naively expect that 
self inductances such that > lO^^^^/Z^ could result in larger conservation probabilities. 

To logarithmic accuracy, L ~ C/c^. From here and Eq. ([6]), tl ~ 
10{e'^klT^ilneW /mc^(^^y/^^j/ D. For the squares in Fig.[6]rL ~ Om^j/D, indicating that 
for quasi ID samples it would be rather difficult to have significantly larger than ^^/-D. 

We now study the evolution of the fiux during a step between t and t + At. We 
write the fiux as {C/2N)J2[^j(t) + ^j(t + ^t)] = Lcl and determine the current by 
NI = J2i^Nj + Isj), where Inj and Isj are the normal and the superconducting current 
at segment j. The normal current is given by I^j = aw[Aj{t) — Aj{t + At) + riAj]/cAt 
and the superconducting current can be taken as Isj = —{2ehwN/mC)lm{ilj*ipj^i), where 
ip{s) = exp[(27ri/$o) Jq ^{s')ds']ilj{s) is the gauge invariant order parameter. From here we 
can obtain 

J2 Mt + ^t) = - At)Aj{t) + 2rLVAj - ^^^^^^Im(^;^,+i) 

(11) 

The fiux enclosed by the loop acts a constraint on the values of for each individual 
segment. It should be noted that the winding number is that of not that of 

Figure E] shows the infiuence of the induced fiux on the conservation probabilities for the 
samples represented by the squares and by the rhombs in Fig. Ej in the range < < 
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FIG. 8: (Color online) Conservation probability ^(-10-^,10-2,6.7^1/1?) as a function of the 
inductive time constant. The squares (black) and the rhombs (blue) correspond to the same 
samples as in Fig. [6l 

1O^0/-D. We found that V does increase with r^, but this increase is moderate and V 
remains far from 1. There is a limited region of steep increase when tl becomes of the order 
of ^"p/D'i the increase in this region is larger for the sample with lower T^. 

VI. INFLUENCE OF THE QUENCHING RATE 



We finally test Postulate 1. For this purpose we varied the temperature with time, 
following the schedule e{t) = e/ tanh(t/rQ). In this way, the sample initially equilibrates at 
e = — ej, quenches through at the rate e = ^j/tq, and finally stabilizes at e = e/. We 
only considered the cases = 10^^ and ej = 10^^; larger values of ej required too short 
iteration steps and smaller values of e/ lead to energy barriers that build up too slowly. 

Our results are shown in Fig. [9l The following conclusions arise from these results: 
(1) Contrary to the expectation from Postulate 1, that (n^) ought to be a function of e/e 
(and therefore independent of ej), {n"^) turned out to be a function of e. (2) (n^) increases 
with e, but for tq < lO^^^^/D this trend saturates and (n^) becomes less sensitive to the 



quenching rate. A similar behavior was observed in Ref. 110|. We may suspect that this 
saturation is related to the decay of the unstable states discussed in Sec. IV A 11 so that 
they don't contribute to the final {n"^)- For samples with weaker superconductivity (smaller 



6/, ^(0) or ^fs), saturation seems to be a stronger effect. (3) For tq > 10 ^l/D, our 
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FIG. 9: (Color online) Variance of the winding number as a function of the quenching rate. The 
initial and final temperatures are Tc(libe/) and the quenching time is tq. Sample parameters, pre- 
sented as (w/lQ-^cm^, 4/10~5cm, r^/K) are □: (5,1,1); o: (10,0.5,0.5); A: (2.5,2,2); *: (1,0.25,5); 
■: (0.25,1,20); -: (40,2,32); •: (0.1,0.5,5.253). The other parameters are as in Sec. HVl except 
for the case of (^y/lO^'^cm^, .£e/10~^cm, Tc/K,e/) = (40,2,32,10^^) in which we took iteration 
time steps of 3.3 x 10^^.^|/D. The black symbols are for e/ = 10^^ and the blue symbols for 
ef = 10^^. The statistical uncertainty of (n^) was estimated as [((f^^) — (n^)^)/(number of real- 
izations for which the evaluation was repeated)] The brown straight lines describe the relation 
(n2) = 9.6 X lQ-\efep/rQDYl^C/{mipY'^- 

results can be described by the expression (n^) = 9.6 x 10"^(e/^^/rQD)^/^C/(^(0)^^)^/^ = 
9.6 X 10~^(e//rQD)^/^C/^(0)^/^ (brown lines in Fig. |9]). In spite of conclusion (1) and the 



differences in the model, we find the scaling (ra^) oc Tq^^^ as in Ref. llOl We note that (n^) 

— 1/2 

is independent of the cross section and proportional to ^ . 



VII. CONCLUSIONS 



We have analyzed the probabilities for final states with permanent currents when a ID 
superconducting loop is quenched through its critical temperature in the absence of magnetic 
flux. The analysis was based on the time-dependent Ginzburg-Landau model with thermal 
fluctuations. 

Our results resemble those predicted by the Kibble-Zurek mechanism (KZM), but mod- 
iflcations are required. The flnal winding number distribution seems to depend on the time 
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derivative of the reduced temperature, e, rather than on the logarithmic derivative e/e. 
For moderate quenching times Tq, the variance of the final winding number is proportional 
to 'Tq^^^ ■ This result is in agreement with KZM, that predicts variances proportional to 
Tq^^^^^''^\ provided that we use the mean field critical exponents z = 2 and v = 1/2. For 
quenching times shorter than 10~^^|/D, the variance of the winding number saturates. Like- 
wise, we found that a large fraction of the "topological defects" that are present immediately 
after the temperature is abruptly lowered below decay during a lapse of time of the order 
of 10~^^^/D. The requirement of a sort of "incubation time" before the winding number 



becomes locked was also noticed in Ref. 



10 



Acknowledgment s 

This research was supported by the Israel Science Foundation, grant No. 249/10. Numeric 
evaluations were performed using computer facilities of the Technion — Israel Institute of 
Technology. I am grateful to Arnab Das, Francesco Tafuri, Vladimir Zhuravlev and Wojciech 
Zurek for their answers to my inquiries. 



Appendix A: Quadratic Approximation for the Distribution of the Order Parameter 



We consider a uniform loop and distinguish between the cases above and below the critical 
temperature. 



1. T > Tc 

In this case an undemanding approximation can be obtained by neglecting and 
CAj/N(^Q in Eq. Decomposing into Fourier components, ipj = V^fcC^'^*-''^/^, we can 

write 

G = J2 (ReVfc + ImVfc) , (Al) 

with 
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From here, the equihbrium distributions for Keifk and Imipk are gaussian, with variance 
kBT/2Ek. 



2. T <Tc 

In this case it will be more realistic to assume that the order parameter is close to the 
ground state, i.e., 

tpj = ^/e/2eiO)P' + uj + ivj , (A3) 

where Uj and Vj are real and small. Substituting Eq. (1A3P into Eq. (jlj) and keeping terms 
up to quadratic in u and v, we are left with 

G = constant + ^^^^(o) + 1^ " + (^i+i " ^i) \ ■ (^4) 

We introduce now the decomposition Uj = Yl!k=iPkG^''^^^^^ = Yl!k=i Ik^^''^^'^^^ with 
PN-k = pI and gAT-A: = ql; in the decomposition for Vj we do not include a constant term, 
since it would just lead to multiplication of the order parameter by a uniform phase. Taking 
N even, the free energy becomes 

G = constant + plfE^ + p^/2^](r/2 + g^/2^]v/2 

N/2~l 

+ J2 [(f^e^p, + Im2p,)E« + {Re\k + l^\k)El] (A5) 



k=l 



with 

From here, the equilibrium distributions for Kepk, Imp^, Reg^ and Img^ are gaussian, with 
variances given by the appropriate ksT /2E^'^ . 
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